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1. Introduction. 

This paper is concerned with a class of time-optimal control problems for the swing 
and the ski. 

We first consider the motion of a man standing on a swing. For simplicity, we neglect 
friction and air resistance and assume that the mass of the swinger is concentrated in his 
baricenter B. Let 6 be the angle formed by the swing and the downward vertical direction 



and r the radius of oscillation, i.e. the distance between the baricenter B and the center of 



rotation O. Assuming that the mass is normalized to a unit, the corresponding Lagrangean 



function takes the form: 
Z ' 1 

L = - (r 2 + r 2 6 2 ) + g r cos(9) (1.1) 

i 

t— i ■ where g is the gravity acceleration. 

> ■ 

We assume that, by bending his knees, the swinger can vary his radius of oscillation. This 

amounts to the addition of a constraint r = u(t), implemented by forces acting on £?, 

parallel to the vector OB. The function u(t) can be regarded here as a control, whose 

values are chosen by the man riding on the swing, within certain physical bounds, say 

"^ ' u(t) G [r_, r_|_] with < r_ < r+. 

Writing (1.1) as a first order system for (6, 9) we obtain an impulsive control system. On 

the other hand, using the coordinates x\ = 9, x 2 = Or 2 (the angular momentum), following 
Ch < 

^ ■ Alberto Bressan (1993), we obtain the nonimpulsive control system: 

c-i— i ■ 

•th : ( ±i - x 2 /v 2 

X I ve[r_,r + }. (1.2) 

[ X2 = —g v sin(xi) 

For the ski, we use the model considered in Aldo Bressan (1991), Bressan and Motta 
(1994), with a special approximation of the skier. 

Our main concern is the existence and the structure of time-optimal controls. As 
a first step, we introduce an auxiliary control system E with control entering linearly 
and establish the relationships between this and the original system. We then perform 
a detailed study of the auxiliary system, using the geometric techniques developed in 
Sussmann (1987 a,b,c) and in Piccoli (1993 a,b). These geometric techniques permit us to 
know the structure of the time optimal trajectories in generic cases and possibly to solve 
explicitely some optimization problem. See section 4 for examples. In turn this provides 
accurate information on the time-optimal controls for the original system (1.2). 
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Since the set of admissible velocities for E is always convex, a time optimal control 
exists for general boundary conditions. However, we prove that a time optimal control for 
(1.2) does exist only when the corresponding optimal control for E is bang-bang. In this 
case the optimal trajectory is the same for the two control systems. We show that, for 
every control constraint v G [r_,r+] with < r_ < r+, there exists some pair of points 
x,x G M 2 such that no control for (1.2) steers x to x in minimum time. 
Then we consider the problem of reaching with the swing a given angle 9 in minimum 
time, having assigned initial condition 6,6 r 2 . This can be stated as a Mayer problem 
for (1.2). We show how to solve this problem under the assumptions \9\, \9\ < 7r/2. To 
solve numerically this problem, it sufficies to find a suitable class of solutions to (1.2) with 
constant control v G {r-,r+}, and then to solve the complementary linear system (2.6) 
for each trajectory. The optimal trajectory is carachterized by the (final) transversality 
condition. We also show that, for a special class of boundary conditions and data r±, no 
optimal solution exists. 

For the ski model we obtain similar results under the assumption that the curvature of the 
ski trail is constant. 
For an introduction to impulsive control systems we refer to Alberto Bressan (1993). 

Acknoledgments. We wish to thank Prof. Alberto Bressan and Prof. Aldo Bressan 
for having suggested the problem and for many useful advises. 



2. Basic definitions. 

By a curve in M n we mean a continuous map 7 : I h^ M n , where I is some real 
interval. We use the symbol Dom for the domain so that if 7 : I 1— >■ M n then Dom('j) = I. 
We use the symbol 7 \ J, where J C Dom('j) is an interval, to denote the restriction of 7 
to J. 

A C 1 vector field on M 2 is a continuosly differentiable map F : M 2 — ► M 2 . It can be 
written in the form: 

F = a d x + /3 d y (2.1) 

where d x , d y are the constant vector fields with components (1,0), (0,1) respectively. If 
we have the representation (2.1) for a vector field F then we use the symbol V-F to denote 
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the 2x2 Jacobian matrix: 

( da da 

dx dy 

d/3 d§_ 

dx dy 

The Lie-bracket of two vector fields F, G is the vector field: 

[F,G] = VG-F - VF-G. (2.2) 

We consider a two dimensional autonomous control system: 

x = h(x,u) uEU (2-3) 

where U C M m is compact and h G C 1 (iR 2 x M m , M 2 ). A control is a measurable function 
u : [a, b] i— ► U where — oo < a < b < +oo. As for the curves we use the symbol Dom for the 
domain. A trajectory for a control u is an absolutely continuous curve 7 : Dom{u) 1— > M 2 
which satisfies the equation: 

7 (t) = fc( 7 (t),u(*)) 

for almost every t e Dom(u). 

If 7 : [a, 6] 1— ► M 2 is a trajectory of (2.3) we use the symbol In(-f) to denote its initial 
point 7(a) and Terrain) to denote its terminal point 7(6). We define: 

T( 7 ) = b - a. 

i.e. T('j) is the time along 7. 

A trajectory 7 is said to be time optimal if for every trajectory 7' with In(Y) = Infa) 
and Term(^') = Term^), we have that ^(7') > T(7). 

If iti : [a, 6] 1— ► [/ and «2 : [b, c] 1— > [/ are controls, we use W2 * «i to denote the control 

defined by: 

(ui(t) for £ G Dom(ui) 
■ 
U2(t) for t G Dom{u2) 

This control is called the concatenation of «i and u^- 

If 71 : [a, b] \— > iR 2 , 72 : [6, c] 1— > iR 2 are trajectories for -ui and u^ such that 71(6) = 72(b), 

then the concatenation 72 * 71 is the trajectory: 

7i(t) for t G _Dom(7i) 

(72*7i)00 

72 (t) for t G Dom{ r )2) 
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Now we consider a control system with control entering linearly: 

x = F(x) +uG(x) \u\ < 1 (2.4) 

where F, G are two C 1 vector field on M 2 and u is scalar. We define: 

X = F-G Y = F + G. 

An X— trajectory is a trajectory corresponding to a constant control u whose value is equal 
to — 1. We define Y — trajectories in a similar way, using the control u = +1 rather than 
u = — 1. An X * Y — trajectory is concatenation of a Y — trajectory and an X— trajectory ( 
the Y — trajectory comes first), and similarly is defined &Y* X— trajectory. 
A bang-bang trajectory is a trajectory that is a concatenation of X— and Y— trajectories. 
A time t G Dom(^) is called a switching time for 7 if, for each s > 0, 7 \ [t — s, t + s] is 
neither an X-trajectory nor a ^-trajectory. If t is a switching time for 7 then we say that 
7(t) is a switching point for 7, or that 7 has a switching at 7(t). 

For the rest of this section we call E the system (2.4). An admissible pair for the 
system E is a couple (u, 7) such that u is a control and 7 is a trajectory corresponding 
to u. We use the symbol Adm(Y) to denote the set of admissible pairs and we say that 
(u, 7) G Adm(E) is optimal if 7 is optimal. 

A variational vector field along (u, 7) G Arfm(E) is a vector-valued absolutely contin- 
uous function v : Dom{^) 1— ► iR 2 that satisfies the equation: 

v{t) = ((VF)( 7 (t)) + «(t)(VG)(7(*))) • «(*) (2.5) 

for almost all £ G Dom(^). 

A variational covector field along (u, 7) G A<im(E) is an absolutely continuous function 
A : Dom(^) i-» iR 2 that satisfies the equation: 

A(t) = -A(t) • ((VF)( 7 (t)) + u(t)(VG)(7(*))) (2-6) 

for almost all t G Dom{^y). Here -R 2 denotes the space of row vectors. In Sussmann (1987 
a) it was proved: 

Lemma 2.1 Let (u, 7) G A<im(E) and Ze£ A : Dom(^) 1— ► iR 2 6e absolutely continuous. 
Then X is a variational covector field along (u, 7) z/ and onfo/ if the function t 1— > A(t) ■ w(t) 
is constant for every variational vector field v along (n, 7) . 
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The Hamiltonian H : Ml x M 2 x M i— ► M is defined as 

H(\,x,u) = \-(F(x)+uG(x)). (2.7) 

If A is a variational covector field along (u, 7) G Adm(E), we say that A is maximizing if: 

H(\(t),i(t),u(t)) = max {H(\(t),i(t),w) : \w\ < 1} (2.8) 

for almost all t G Dom{^f). 

The Pontryagin Maximum Principle (PMP) states that, if (u, 7) G Aim(E) is time 
optimal, then there exists: 

(PMP1) A non trivial maximizing variational covector field A along (u, 7) 

(PMP2) 4 constant A < suc/i tfm£ : n(X(t),^y(t), u(t)) + A = for almost all £ G Dom^). 

In this case A is called an adjoint covector field along (u, 7) or simply an adjoint 
variable, and we say that (7, A) satisfies the PMP. 

If A is an adjoint covector field along (u, 7) G Adm(Y,), the corresponding switching 
function is defined as: 

A (t) = A(t)-G( 7 (t)). (2.9) 

From the above definition it follows: 

Lemma 2.2 Let (u, 7) G Adm(E) 6e optimal and let A 6e an adjoint covector field along 
(u, 7) . Then: 

a) The switching function <f>\ is continuous. 

b) If (pit) > for all t in some interval I , then u(t) = 1 for almost all t G I and 7 \ I is 
a Y -trajectory. 

c) If <pit) < for all t in some interval I , then -u(t) = —1 for almost all t G I and 7 |" I 
is an X -trajectory. 



For each x G M 2 , one can form the 2x2 matrices whose columns are the vectors F, G 
or [F, G\. As in Sussmann (1987 a), we shall use the following scalar functions on IR 



1 **-" ? 

2. 



A A (x) = det(F(x), G(x)) (2.10) 

A B (x) = det(G(x), [F, G]ix)) (2.11) 



where det stands for determinant. 



6 B. PICCOLI 

Consider (u,j) G Adm(T,), to G Dom^) and vq G iR 2 . We write v(vQ,to;t) to denote the 
value at time t of the variational vector field along (u, 7) satisfying (2.5) together with the 
boundary condition v(to) = vq. In Piccoli (1993 a) it was proved the following: 

Lemma 2.3. Let (u, 7) G Aim(E), to G Domi^y), and vq G iR 2 , vq 7^ 0. For every t snc/i 
t/iat G(j(t)) j£ 0, define the angle: 



a(t) = argU, u(G( 7 (t)),t; t )J, (2.12) 

Then, one has: 

sgn(d(t))=sgn(A B ( 7 (t))). (2.13) 

A point x G -R 2 is called an ordinary point if 

A^OrO-AsOr^O. 

On the set of ordinary points we define the scalar functions /, as the coefficients of the 
linear combination 

[F, G](x) = f(x)F(x) + g(x)G(x). (2.14) 



By direct calculations we have: 



/= " B 



A point x at which Aa(x)Ab(x) = is called a nonordinary point. A nonordinary 
arc is a C 1 one-dimensional connected embedded submanifold S of iR 2 , with the property 
that every x G S is a nonordinary point. A nonordinary arc will be said isolated, and will 
be called an INOA, if there exists a set Q satisfying the following conditions: 

(CI) O is an open connected subset of iR 2 

(C2) S is a relatively closed subset of O 

(C3) If x G O — 5 then x is an ordinary point 

(C4) The set Q — S has exactly two connected components. 

A turnpike is an isolated nonordinary arc that satisfies the following conditions: 

(51) For each x G S the vectors X(x) and Y(x) are not tangent to S and point to opposite 
sides of S 

(52) For each x G S one has Ab(i) = and Aa(x) 7^ 
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(S3) Let O be an open set which satisfies (C1)-(C4) above. If Ox and fiy are the connected 
components of O — S labelled in such a way that X(x) points into Ox and Y(x) points 
into fly, then the function / in (2.14) satisfies 

f(x) > on Q Y f(x) < on Ox- 

Next, consider a turnpike S and a point xq G 5". We wish to construct a trajectory 
7 of (2.4) such that 7(^0) = xq and 7(t) G S for each t G Dom(^) = [to,ti\. Clearly, one 
should have As (7(2)) = for all t. Since A B (7(^0)) = 0, it suffices to verify that: 

j t A B { 7 (t)) = (VA B ^)(t) = 0. 
The above holds provided that 

(VA B ■ uG)(<y(t)) + (VA B • F)( 7 (t)) = 0. 

Assuming that 

(VA B • G)(x) ^0 VxeS, (2.15) 

the values of the control u are thus uniquely determined by 

. VA B ■ F(x) 
" = #*)= - VA B .G(xY (2 - 16) 

A turnpike is regular if for every x E S (2.15) holds true and \<f>(x)\ < 1. A trajectory 
7 is said to be a Z-trajectory if there exists a regular turnpike S such that {^(t) : t G 
Dom{^)} C S. 

An isolated nonordinary arc or INOA S is a barrier in O (where O satisfies (C1)-(C4)) 
if it verifies: 

(SI') For every x E S, X(x) and Y(x) point to the same side of S 

(S2') Each of the function A^, A B is either identically zero on £ or nowhere zero on S. 

A point x G M 2 is a near ordinary point if it is an ordinary point or belongs to an INOA 
that is either a turnpike or a barrier. 

Remark 2.1. 

The definition of near ordinary point given in Sussmann (1987 a) is more general, but 
for our pourpouses this simpler definition is sufficiently general. 



8 B. PICCOLI 

3. Preliminary Theorems. 

In this section we show three theorems on control systems with control appearing 
linearly as in (2.4) and prove a theorem relating the control systems for the swing and the 
ski with this type of control systems. 

Consider the control system (2.4). For ordinary points we have the following: 

Theorem 3.1 Let O C M 2 be an open set such that each x G O is an ordinary point. Then 
all time optimal trajectories 7 for the restriction of (2.4) to O are bang-bang with at most 
one switching. Moreover if f > throughout Q then 7 is an X,Y or Y * X — trajectory, if 
f < throughout Q, then 7 is an X,Y or X * Y— trajectory. 

For the proof see Sussmann (1987 a, p. 443). For near ordinary point we have a similar 
theorem on the local structure of optimal trajectories, see Sussmann (1987 a, p. 459): 

Theorem 3.2 Let x be a near ordinary point. Then there exists a neighborhood O of x 
such that every time optimal trajectory 7 for the restriction of (2.4) to O is concatenation 
of at most five trajectories each of which is an X—, Y — or Z— trajectory. 

Let consider now a control system as in (2.3): 

x = h(x,u). (3-1) 

If W\,W2 G -R 2 , define the triangle: 

C(w\, W2)={w G IB 2 : w = Xwi + [iu)2\ A, /j, > 0; A + \x < 1} 

and consider the condition: 

(PI) There exist two C 1 vector fields w ± {x) that for every x either are linearly independent 
or have the same versus, and such that: 

{w ± (x)} C {h(x,u) :ueU} 

{h(x,u) : u G U} C {x : x = Xw + (x) +fiw~(x). A, \x > 0, < A + /U < l}U{w (x)} C 

C C(w + (x), w~(x)). 
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Suppose that (3.1) verifies (PI). Then we can define a system of the form (2.4) choosing: 



f = ^±m, g = ^m (3.2) 

thus (F ± G)(x) = w ± (x). We can also define the map: 

V : M 2 x U -»■ [-1, 1] (3.3) 

in the following way. If w ± (x) are independent then V(x, v) — u if and only if h(x, v) and 
F(x) + uG(x) are parallel. Otherwise, V(x, v) is constantly equal to +1 if F(x), G(x) have 
the same versus (that is if \w + \ is bigger than \w~\) and to —1 if the opposite happens. 
From (PI) we have that this map is well defined. 

Given two points x,x G M 2 we consider the two endpoints problem of steering x to x in 
minimum time. Let define the value functions: 

V(x,x) = inf {T(-f) : 7 trajectory of (2.4), (3.2), 771(7) = x,Term(^) = x} 

V(x,x) = inf {T(7) : 7 trajectory of (3.1), 171(7) = x,Term(^) = x} 
We have the following: 

Theorem 3.3 Assume that every point is near ordinary for (2.4), (3.2) and consider 
two points x,x G IB? : 

(i) there exists a trajectory 7 of (2. 4), (3. 2) such that T{^) = V(x,x) 
(ii) if Y is the set of trajectories of (2. 4), (3. 2) corresponding to bang-bang controls and 

V'(x, x) = inf{T(-f) : 7 G T, In(-y) = x, Term{p/) = x} then: 

V'(x,x) = V(x,x). 



Proof. Assume that V(x,x) < +00, otherwise there is nothing to prove. The state- 
ment (i) follows from the convexity of the set of velocities for (2. 4), (3. 2), see for example 
Lee and Markus (1967). 

Let 7 be a time optimal trajectory steering x to x. From Theorem 3.1 and 3.2 we have that 
7 is concatenation of X, Y and Z— trajectories. If 7 is bang-bang we are done. Assume 
now that 7 \ [to, £1] is a Z— trajectory. From the definition of turnpike we have that for 
every t G [to, £1] the two vectors X(7(t)), Y(-f(t)) points to opposite sides of the image of 
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7. If 7(f) is sufficiently near to 7(to) then we can construct a bang-bang trajectory 7 with 
one switching steering 7(^0) to 7(^1), in the following way. For \t — to\ sufficiently small the 
Y— trajectory 7 + passing through 7(^0) an d the X— trajectory 7" passing through 7(f) 
meet each other in at least one point. Let x be the first point in which 7+ intersect the 
image of 7", after passing through 7(^0) • We can construct 7 following 7+ up to the point 
x and then 7" up to the point 7(t). Divide now [to, ii] into n equal subintervals inserting 
the points ki=to + (i/n){ti — to), i = 0, . . . ,n. If n is sufficiently large we can construct, 
as above, a bang-bang trajectory steering 7(fcj) to 7(^+1), i = 0, . . . , n — 1. Let 7 n be the 
concatenation of these trajectories. Let K be a compact neighborhood of the image of 7 
and let M=2 maxK {\F\ + \G\) < +00. If n is sufficiently large then by construction 7^ 
lies in K, hence: 

n n 

The subset J of Domini) on which 7 is a Z— trajectory is a finite union of closed interval. 
Otherwise we can constuct a sequence t n in d J converging to a time t G Dom( , y) and then 
7(t) is not a near ordinary point. Repeating the same reasonings for every subinterval of 
J we obtain a sequence of bang-bang trajectories converging uniformly to 7. 

Notice that every bang-bang trajectory for (2. 4), (3. 2) is a trajectory for (3.1), then 
V(x,x) < V'(x,x). We can prove the following: 

Theorem 3.4 Assume that (PI) holds true and that every point is a near ordinary point 
for (2. 4), (3. 2). Consider two points x,x G M 2 then 
(i) if there exists a bang-bang time optimal control u for (2. 4), (3. 2) steering x to x, then 
there exists a time optimal control v for (3.1) corresponding to the same trajectory 7 
ofu, i.e. h( 1 (t),v(t)) = F( 1 (t)) + u(t)G( 1 (t)) G {^(7^))}; 
(ii) if every time optimal control u for (2. 4), (3. 2) is not bang-bang (i.e. if j(t) belongs 
to a turnpike for some interval) then the time optimal control for (3.1) does not exist 
but we have: 

V(x,x) = V(x, x) 

and for eache > there exists a control v, corresponding to a trajectory n steering x to 
x, such that h(i](t), v(t)) G {w ± (v(t))} f or eac h t *= Dom{^) and T(rj) < V(x, x) + e. 

Proof. Suppose first that there exists a bang-bang time optimal control u and, by 
contradiction, that there exists v control of (3.1), with corresponding trajectory r\ steering 
x to x, such that T(tj) < T("f). If h(i](t),v(t)) G {w ± (rj(t))} for almost every t then r\ is 
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a trajectory of (2. 4), (3. 2) contradicting the optimality of 7. If this is not true let define 
the feedback control u(rj(t))=V(r}(t), v(t)), where V is the map in (3.3). There exists a 
trajectory 7 corresponding to u that is a reparametrization of rj and is a trajectory of 
(2.4), (3.2). Finally \F(r]{t)) + u(r](t))G(r]{t))\ > \h{r]{t),v(t))\, therefore: 

T( V ) > T(7) > T( 7 ) 

that gives a contradiction. 

Suppose now that every time optimal control u is not bang-bang. From Theorems 3.1 
and 3.2 we have that there exists / C Dom^) such that 7(1) is contained in a turnpike. 
Consider a control v for (3.1) corresponding to a trajectory r\ steering x to x. Let define u, 7 
as above and define S = {t : \h(j]{t),v(t))\ < \F(rj(t)) + u(r](t))G(r](t))\}. If meas(S) = 
then 77 is a bang-bang trajectory of (2.4), (3.2) and is not time optimal. If meas(S) > 
then there exists n and 5 > such that 

meastlt: \h(r}(t),v(t))\ < \F(rj(t)) + u(r}(t))G(r}(t))\ - -1 J > 5 > 0. 

Therefore: 

T( V ) > T(7) + 5 - > T( 7 ) 
n 

for every time optimal 7. 

From Theorem 3.3 it follows the second part of the statement (ii). 

Theorem 3.4 shows the relationships between the two control system (2. 4), (3. 2) and 
(3.1). If we are able to determine the time optimal control u for a given problem for 
(2. 4), (3. 2) then we immediately know all about the same problem for (3.1). Indeed if the 
time optimal control u is bang-bang then the trajectory 7 of u is a trajectory also for (3.1), 
is optimal and corresponds to a control v taking values in {r_, r+j. If every time optimal 
control is not bang-bang then the time optimal control for (3.1) does not exists. 
Thank to Theorems 3.1 and 3.2 in most cases we are able to know the local structure of time 
optimal control for (2. 4), (3. 2) and then for (3.1). Using the same methods of Sussmann 
(1987 a,b,c) and Piccoli (1993 a,b), we are able to solve explicitely many optimization 
problems for the swing and the ski models. Some examples will be given in the following 
sections. 



4. Time Optimal Control of the Swing. 



12 B. PICCOLI 

In this section we treat the problem of time optimal control for the swing considering 
the minimum time problem and a Mayer type problem. 

Recall equation (1.2). We have to verify that this control system satisfies the condition 
(PI) of section 3. If sin(xi) = then xi = for every v G [r_, r+] and the set of velocities 
is a segment. If sin(xi) ^ 0, from the second equation of (1.2) we obtain: 

X2 (4.1) 



g sin(xi) 
and replacing (4.1) in the first equation of (1.2): 



Xl = x *9 28 *f(*i). (4.2) 

x 2 

It is easy to check from (4.2) that the set of velocities for (1.2) at a given point lies on a 
branch of hyperbola. Defining w = h(x,r^), (PI) holds. Notice that with this definition 
we have: 

(F + G)(x) =h(x,r_) (F-G)(x) = h(x,r + ). 

Therefore if u is a bang-bang control for the auxiliary system (2. 4), (3. 2) with corresponding 
trajectory 7, then 7 is also a trajectory for (3.1) and corresponds to the control: 

v(t) = r_ if u(t) = 1, v(t) = r + if u(t) = -1. 

Hence we can compute the linear system (3.2) associated to this control system as in 
section 3, explicitely: 

r%+r 2 _ \ / r 2 + -r 2 __ 

2r 2 + ri X2 I G= I 2r 2 + ri X<2 

-| (r+ + r_) sin(xi) J \ | (r+ — r_) sin(xi) 

For simplicity we define: 

a = (r+ r_) 2 b = r + + r_ c = r+—r- d = r\ + r_ (4.3) 

then: 

F = ( 2^ X2 \ G — ( 2^ X2 ] (4.4) 

\—^bsin(xi) J \ § c sin{x\) J' 

To investigate the local structure of time optimal trajectories we have to compute the 
functions Aa,Ab, defined in (2.10,11): 

. gc(d + b 2 ^ 



Aa 



X2 sin(xi) (4-5) 
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[F,G\= ( ~ 9C \* b2) gm(xi) 

\ 9C 4^ b X2 cos ( x l) 

hence: 

Ab = I - x\ cos(xi) + g sin 2 (xi) I . (4-6) 

8 a \ a " / 

Every turnpike is subset of A^ 1 (0), then we have to solve the equation Ag(x) = that is 
equivalent to: 

- x 2 cos(xi) + g sin (xi) = 
a 



that gives: 



''_• 



g a sin 2 (xi) 



b cos(xi) 

By periodicity we can restrict ourselves to the case x\ G [0, 27r]. There is the isolated 
solution (0, 0) and if x\ G]^, ^-[ we have the two solutions: 



, g a sin 2 (x 1 ) 

otherwise there is no solution. 

The two branches of solutions form two curves, one contained in the first quadrant, the 
other in the forth quadrant. The two curves meet each other at the point (it, 0) and 
(0, 0), (tt, 0) are the only points of A" 1 ^) that verify VA B (x) = 0. For each x G A" 1 ^) \ 
{(0, 0), (tt, 0)} we can compute the control <f>(x) defined in (2.8): 

, , . , . . (2 b 2 - d) cos 2 (xi) - d ,, , 

b c (3 cos z (xi) + 1) 

From (4.8) we have: 

rp^ _l_ rp^ 

< -1 



lim 0(xi) = 

1_ > "2 


lim <f>(xi) ■ 


r\ + r 2 _ 


rp^- rp^ 


lim 0(xi) 

X\ — >7T 


2b 2 -2d 
Abe 


r+ r_ 
— ^ > 

rp£ rp^> 


d(f)(xi) 


4 (b 2 + d) 


sin(xi) cos(xi) 



dx\ b c (3 cos 2 (xi) + l) 2 

then <p is increasing for x\ e\tt/2, ix[ and decresing for x\ E]n, 3 7r/2[. Therefore there exist 
£i > 0,£ 2 > such that \(j>(xi)\ < 1 for x x g]tt/2 + e 1? n - e 2 [ U ]7r + e 2 ,3 tt/2 -£i[# 0. 
We have that e 2 = if and only if: 



r+ — r 



=r<! 
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i.e. if and only if: 

^ 1 + V5 

+ ~ 2 

Hence regular turnpikes do exist. We can choose two points x,x G M 2 such that the only 
time optimal trajectory for (2.4), (4.4) that steers x to x is not bang-bang: it is sufficient 
to take two points of the same turnpike, see Sussmann (1987 a). In this case the time 
optimal control for (1.2) does not exist and the second part of Theorem 3.4 applies. 

Now, Theorems 3.1 and 3.2 determine the local structure of the time optimal trajec- 
tories for the system (2. 4), (4. 4) and then for the swing model. 

Observe that A^ 1 (0) = {(xi,^) : ^2 = or sin(xi) = 0}. Therefore it easy to check 
that A^ (0)\{(0, 0), (tt, 0)} is union of a finite number of INOAs that are barriers. More- 
over, from the reasoning above we have that Ag (0)\{(0, 0), (tt, 0)} is union of a finite 
number of INOAs each of which is either a turnpike or a barrier. Hence every point of 
O=iR 2 \{(0, 0), (7r,0)} is near ordinary. Notice that if 7 is a trajectory of (2. 4), (4. 4) and 
7(£) G {(0,0), (-7T, 0)} for some time t G Dom^), then 7 is a constant trajectory. Indeed, 
in this case, the control has no effect being X(-f(t)) = Y(j(t)) = 0. Moreover, (0, 0), (ir, 0) 
are the only point in which either X or Y vanishes. Thus if 7 is a trajectory of (2. 4), (4. 4) 
we have that either 7 is constant or ^(t) G O for every t G Dom{^). In the latter case we 
can apply Theorems 3.1 and 3.2. 
Given the expression of Aa, Ag we can calculate / of (2.14): 

. c \^x\ cos(xi) +g sin 2 (x] 
f{x 1 ,x 2 ) = — ^ 



2x2 sin(xi) 

Consider the set: 

Q=|(xi,x 2 ) : \xi\ < -| 

and a time optimal trajectory 7 that verifies 7(t) G Q for every t G Dom{^). We have, 
from Theorems 3.1 and 3.2, that 7 is bang-bang and that 7 can switches from control 
+1 to control —1 if x<i sin(xi) > and from control —1 to control +1 if x<i sin(xi) < 0. 
Therefore Q\A~± (0) is divided is four parts and on each parts only one kind of switching 
is permitted. This correspond to the fact that if the swing is raising his distance from the 
earth then the swinger can change only from control r_ to control r+, instead if the swing 
is lowering his distance from the earth then the swinger can change only from r+ to r_. 
In this case, the map V of (3.3) is bijective for almost every x and we can establish a 
bijective correspondence between the trajectories of (1.2) and of (2. 4), (4. 4). In particular 
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if 7 : [a, b] — ► M 2 is a trajectory of (2. 4), (4. 4) then there exists a trajectory n : [c, d] — > iR 2 
of (1.2) verifying 77(c) = 7(a), rj(d) = 7(6) and ?y(t) G 7([o, &]) for every t G [c, d]. Therefore 
also some geometric properties of (2. 4), (4. 4) hold for (1.2). For example the reachable sets 
from a given point are the same for (2. 4), (4. 4) and for (1.2). 

In Alberto Bressan (1993) it was considered the problem fo raising the amplitude of the 
first half oscillation starting from a given point x G M 2 . Notice that if sin(xi) > then 
F 2 — G 2 < F 2 + G 2 < and the opposite happens if sin(a;i) < 0. Hence, comparing the 
vector X, Y, it is easily seen that if we want to reach the maximum amplitude of the first 
half oscillation we must choose the control +1 if x 2 sin(xi) > and —1 if x 2 sin(xi) < 0. 
For the swing this means to choose r_ if we are raising the distance from the earth and r + 
in the other case. Obviously this is also the optimal control to raise the amplitude after a 
given number of oscillation. 

Using the Pontryagin Maximum Principle and Lemma 2.1,2.2 and 2.3 we can prove 
the following: 

Lemma 4.1. Assume that r\ : [a, b] — ► M 2 is a bang-bang trajectory of (1.2), rj(t) G Q for 
every t G [a, b], and that there exists t\,t 2 G (a, b) such that either the first or the second 
component of rj(ti) vanishes and the same hold for n(t 2 ). If rj has no switching then r\ can 
not be optimal. 

Proof. Assume for example that rj corresponds to the constant control +1 and, by 
contradiction, that rj is optimal. Observe that r\ is a trajectory of (2. 4), (4. 4), then if n 
is optimal there exists an adjoint covector field A(t), t G [a, b], along 77. We have, from 
Lemma 2.2, that A(rji) • G (r}(h)) > and, from Aa(t/(*i)) = 0, that G(7/(ti)),y(?/(ti)) are 
parallel. From Lemma 2.3 and (4.6) the function: 



ft 



t(*)=ar(/(y(77(ti)), u(G(7/(t)),t;ti)) 

is strictly increasing. It is easy to check from (2.5) and (4.4) that the function: 

ip(t)=det(v(G(r](t2)),t2;t), v{Y{r ] {t 2 )),t 2 -t) 



is constant. Assume for example that t 2 > t\. We have that G(t 2 ) = v(G(rj(t 2 )),t 2 ; t 2 ) and 
Y(t 2 ) = v(Y(r](t 2 )),t 2 ;t 2 ) are parallel because A a (vi^)) = 0. Hence v(G(r](t 2 )),t 2 ;ti) 
and v(Y(r}(t 2 )),t 2 ',ti) are parallel, but from (2.5) it follows v(Y(rj(t 2 )),t 2 ;ti) = Y(rj(ti)). 
We conclude that a(t 2 ) = kir for some integer k > 0. 
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Now, from Lemma 2.1 and 2.2: 

X(t) ■ G{r}{t)) = \(t 1 )-v(G( V (t)),t;t 1 ) >0 VtG [a, b]. 

But since a(t 2 ) > ^ the vector v(G(r}(t)),t; tA for t G [a, b] makes a rotation of an angle 
strictly greater than 7r and then there is no vector A(£i) for which the above inequality can 
hold. 

It is useful the following: 

Lemma 4.2. Assume that 7 is a time optimal trajectory of (2.4), (4-4)? that 7 has a 
switching at time t\ G Dom{pi) and that A ^4 (7(^1)) = 0. Then A ,4 (7(^2)) = 0, t 2 G 
Domi^j), if an d on ^V it '^2 is a switching time for 7. 

Proof. Let u(t), t G Dom(^), be the control corresponding to 7. Since 7 is optimal 
there exists an adjoint covector field A along (u, 7) . Let £2 be the first time that either is 
a switching time or that verifies A ^4 (7(^2)) = 0. We have that u \ \t\,t 2 \ is constant, say 
u = 1, and that A(*i) • G(7(*i)) = 0. But Gfaftfci)) and ^(7(^1)) + G(7(ti)) are parallel, 
hence A(ti) • [^(7(^1)) + G(7(ti))] = 0. From Lemma 2.1 have: 

A(t 2 )- [F( 7 (t 2 ))+G( 7 (t2))] = X(t 1 )-v(F( 1 (t 2 )) + G( 1 (t 2 )),t 2 ;t 1 ) = 

= A(t 1 ).[F( 7 (t 1 )) + G( 7 (ti))]=0. 

Now if A ,4 (7(^2)) = we have that G^^X-F^fe)) + ^(7(^2)) are parallel then t 2 is 
a switching time. On the other hand, if t 2 is a switching time then A(t2) • ^(7(^2)) — 0, 
hence being X(t 2 ) ^ we have that ^(7(^2)) and ^(7(^2)) + G("f(t 2 )) are parallel. This 
means A^i^fe)) = 0. We can argue in the same way for the other switching times. 

Consider now the problem of reaching, with the swing, an angle 9 G (0, 7r/2) in 
minimum time with given initial condition. We argue about the corresponding problem 
for (2. 4), (4. 4) with initial point x G Q. We can solve this problem constructing all the 
time optimal trajectories starting from x and using the (final) transversality condition of 
the PMP, see for example Lee and Markus (1967), to select among these trajectories the 
optimal one. To construct all the time optimal trajectories we can proceed as in Piccoli 
(1993 a,b). 

Assume, for example, that x = (27, £2), %\ > 0, x 2 > 0, being similar the other cases. 
Let 7 be a time optimal trajectory that verifies In('j) = 7(0) — x. If 7 lies on the first 
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quadrant for t G [0, ti] then 7 \ [0, t{\ is a Y — trajectory or a X * Y— trajectory. Let 7+ 
be the Y — trajectory satisfying 7 + (0) = x, and let t + be the first positive time, if any, for 
which 7 + (t + ) belongs to the x\— axis. 

First we want to prove that if t + < 00, 7(0) = x and 7 \ [0, t + + e] is a 1"— trajectory for 
some e > 0, then 7 is not time optimal after a given time. Assume, by contradiction, that 
7 is time optimal. From Theorem 3.1, being / > in the fourth quadrant, we have that 
necessarily 7 is a Y — trajectory until it reaches the £2— axis, say at time t\. Now 7 can 
not switch at time t\ because of Lemma 4.2. Hence we can apply Lemma 4.1 obtaining a 
contradiction. 

Then if we want to contruct all the time optimal trajectories it is enough to consider the 
trajectories 7 S , 7 S (0) = x, that follow 7 + for a given time s G [0,t + ] and then switch 
to control — 1. Now if 7 S is time optimal then there exists an adjoint covector field A s 
along 7 S . Since X s (s) ■ G("f s (s)) = we can determine A s up to the product by a positive 
scalar. Hence we can compute the switching function <p s = X s ■ G7(7 S )) and, by Lemma 
2.2, determine the behaviour of 7$, that is its switching times. Using Lemma 4.2 we can 
see that, for every s, 7 S has to switch on each quadrant. More precisely 7 S will make a 
second switching before reaching the X2— axis, then a third switching before the second 
time of intersection with the x\— axis and so on. The set of switching points of 7 S form 
some manifolds called switching curves, see Piccoli (1993 b) for the exact definition. Each 
switching curve lyes on a quadrant. After a given number of oscillations some 7 S reach the 
manifold M.={{x\,X2) '■ x\ = 6}. To select the optimal trajectory between these 7 S , we 
can use the transversality condition that, in this case, is X s (t s ) • (0, 1) = if t s is the first 
time of intersection of 7 S with M.. 

To solve numerically our problem, first we have to solve the equation for the swing (1.2), 
for constant control r+ and r_, and a class of initial data (that is to approximate some 
elliptic integrals) . Indeed every 7 S is a finite concatenation of such trajectories. Then we 
have to consider the complementary equation (2.6) for j 8 . This is a linear system and can 
be solved numerically by usual methods, determining, up to a scalar, the adjoint covector 
field A s with initial condition X s (s) ■ G( , y s (s)) = 0. Finally the transversality condition 
determines a value s, and then the corresponding trajectory 75. 

If x (£ Q and \9\ > tt/2 then the construction of time optimal trajectories is more difficult. 
We can do it following Piccoli (1993 a,b), but in this case we have to take into account the 
turnpikes. It can happen that some time optimal trajectories are Z— trajectories for some 
time interval of positive measure, and the optimal control for the swing does not exist. An 
example of this situation is given in the following. 



18 B. PICCOLI 



We now consider the Mayer problem: 

f-y(t) = F( 7 (t)) + «(t)G( 7 (t)) 

^ 7(0) = x (4.9) 

[ min {£ : 3 7 such that 7 (t) G 7W} 

where M.={{x\, £2) :xi = c},cG iR. We consider the case in which c = 3 7r/2, x = (xi, X2) 
belongs to a turnpike and x\ e]tt, Stt/2[, x^ > 0. Our aim is to show the existence of some 
values of r± such that the problem (4.9) with the dynamics (1.2) does not have a solution. 

If 7 is a solution then in particular 7 is time optimal for steering x to Term (7) and 
hence is a concatenation of X—, Y — and Z— trajectories. Following the algorithm in Piccoli 
(1993 a,b), we can cover a region of the plane with the time optimal trajectories starting 
from x. Let S be the turnpike to which x belongs. The region of first quadrant below 
S is covered by Y — trajectories originating from S and the region over S is covered by 
X— trajectories. If x' is the endpoint of S that comes after x for the orientation given 
by X, Y, then the X— trajectories that cross A^ (0) over x' must switch changing the 
control to +1. The points in which these trajectories change control form a curve called 
switching curve, see Piccoli (1993 b). Therefore if 7 is a time optimal trajectory we have 
two possibilities: either 7 = 72 * 71 where 71 is a ^-trajectory possibly trivial and 72 is 
a Y— trajectory, or 7 = 73 * 72 * 71 where 71 is a Z-trajectory possibly trivial, 72 is an 
X— trajectory and 73 is a Y — trajectory. 

Our aim is to find some r± such that the solution to (4.9) is of the first type with 71 
not trivial, i.e. Dom('j) is not a single point. 

Let define 7+ as the Y— trajectory that verifies 7n(7+) = 7^(0) = x and let 7^ be 
the Z— trajectory that satisfies i"n( 7 f ) = 7 f (0) = x. We have that 7+ satisfies: 

xi = ^fx 2 



hence: 



We define: 



and u> by: 



g (c—b) ■ / \ 

X2 = ^2 — sin[x\) 



^(^)^f!^)fiS^ )Mxi) . (,10) 



a = ^±^, p= 9(c-b) (4n) 

2a ' 2 v ' 



oo 2 = - a p. (4.12) 
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We can solve (4.10) in t using the first integral: 

x\ — 2 ui 2 cos(xi) 



and obtaining: 



•2 2 —2 a 2 / • 2 1 • 2 1 

x t = a x 2 + 4 uj [ sin — sin — 



Now let x=7^~(T) where T is the first time at which 7^" intersects A4 and define: 

4uj 2 



k = sin — 



a 2 x\ + 4 uj 2 k 2 ' 



After straightforward calculations from (4.13,14) we obtain: 

1 r 1 dx x 



T = 



V 7 " 2 x 2 2 + 4uj 2 k 2 J Sl 1 1 _ a 2 sin 2 x± 



and, using the substitution 0=^: 



d,e 



'a 2 x\ + Au 2 k 2 J^± a/1 - a 2 sin 2 6 
2 



'a 2 xl + Aui 2 k 2 



Etao,^)-E(ao,^ 



where: 



E(i,e) 



do 



(4.13) 



(4.14) 



(4.15) 



(4.16) 



/o Vl - I 2 sin 2 9 

is the elliptic integral of the second type. Now we want to compute the time along the 
trajectory "f s that satisfies 7 S \ [0, s] = 7^ |" [0, s] and that is a Y— trajectory after the time 
s. Define T(s) to be the first time in which 7 S intersects M.. We have: 



T{s) 



s + 



where: 



a 2 x 2 (s) + Auj 2 k 2 

,- s / \ 7 - - / X 1 [ S 1 

x{s) = 7j(sJ k s = sin 



E a 



3tt 



Ela 



xi(s) 



(4.17) 



4uj z 



2 J s a 2 x 2 (s)+4uJ 2 k 2 ' 
For any s > we can calculate the difference between T and T{s). It is clear that if: 



(4.18) 






(4.19) 



s=0 
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is negative then 7^" is not optimal for (4.9). In computing (4.19) we will use the approxi- 
mations: 

xi(s) — xi + (Fi + 4>Gi)(x) s + o(s) x 2 (s) = x~2 + (F2 + (j)G2){x) s + o(s). 



After some calculations from (4.15,16,17,18) it follows: 






s=0 



l+(a 2 xl + Au 2 k 2 ) 2 A 



■•■ 1 
2 



(l - al sin 2 9) 2 + (l - al sin 2 9) 2 a 2 sin 2 9 d9 



+ 



d + 4>(xi) b c 
2^ 



x 2 [(a 2 x 2 . + 4 u 2 k 2 ) (l-al k 2 )] 



(4.20) 



where: 



2 /,/- n ,x - ■ - , 2- ■ x i x i (d + (j)(xi)bc 

A = a g {<p{xi) c — 0) X2 sin x\ + 2 u> X2 sin — cos — 



age 
2a 



(d + b 2 ) (4>(xi) — 1) X2 sin x\. 



(4.21) 



Notice that the second and the third terms on the righthand side of (4.20) are negative. 
If r_ tends to zero then the second term tends to: 



cos(xi) (3cos (xi) — cos(xi) + 2) (x\—x\ 



sin(xi) 



3cos 2 (xi) + 1 



(4.22) 



Now if we choose r+ sufficiently large then |0(xi)| < 1 for x\ e]tt, tv + e[ and some e > 0, 
hence we can take x\ arbitrarily near to ix. But as x\ tends to n, the expression in (4.22) 
tends to minus infinity. In particular we can choose r_, r_|_ and x\ in such a way that: 



ds 



<0. 



(4.23) 



s=0 



Now it is easy to see that, choosing r±, X\ as above, the Y * X— trajectory leaving from x 
can not be optimal for (4.9), in fact X(x) tends to zero. On the other hand for (4.23) the 
Y — trajectory starting from x can not be optimal. Hence the time optimal trajectory for 
(4.9) contains a non trivial Z— trajectory. Using Theorem 3.4 we obtain that the problem 
(4.9) with the dynamics (1.2) has no optimal solution. 
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Given the time optimal control u for (4.9), corresponding to the trajectory 7, we can 
consider the feedback control v such that, see (3.3): 

P(7(t),v(7(*))) = «(*)■ 

Indeed, in this case, the function V, defined in the third section, is bijective outside the 
coordinate axes and the line {x\ = n}. We have: 



»«'» = «-;$& < 4 ' 24 > 

The trajectory r\ corresponding to the control v is a reparametrization of 7. 

Since 7i(t),?)i(t) > for every t we can consider the inverse functions t^(xi),t v (xi) and 

define 7(^1) = 7(t 7 (^i)), r}{x\) = r]{t fl {xi)). If x = Term(^) then: 

r(7) = r d n 

J Sl F x (7(^1)) + u(t 7 (xt)) d (7(^1)) 
= r 1 v 2 (r ] (x 1 ))dx 1 

Then we can calculate the difference of times T{pf) — T{rj) computing the difference of the 
two integrals. We can obatin a better performance using a bang-bang control but it is 
not easy to compute explicitly one such control and its total variation tends to infinity as 
its time tends to the minimum. The control v has the advantage of being defined by the 
explicit formula (4.24) and of having a fixed variation. 



5. Time Optimal Control for the Ski. 

In this section we deal with time optimal control problems for the ski. We consider a 
skier on a one-dimensional trail described by a curve (x(s),y(s)) in the plane. The skier 
can choose the height u of his baricentre, that is the distance from the point of contact with 
the trail, within certain physical bounds. We refer to Aldo Bressan (1991) for assumptions 
and notations. The dynamics is given by the system: 

s = p 

X{s,u) 

u G [r_,r + ] (5.1) 

P = y^P 2 -M g{l-cu)y'{s) 
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where s is the arclength of the trail, p its conjugate momentum, c(s) the curvature of the 
trail, I is c 2 times the inertial moment of the pair ski-skier with respect to the centre of 
curvature, X s its partial derivative with respect to s, g the gravity acceleration and M is 
the total mass of the pair ski-skier. We approximate the body of the skier with a system 
of n equal masses with height i h/n, i = 1, . . . , n, where h is the height of the top of the 

skier. We have: 

m + m s n 
h = bu, 6 = 2 



"? 



m n + 1 



and, see Bressan (1991): 

n / . v 2 

X(s, u) = ac 2 (s)+m s ^ V ( 1 - c(s) bu - (5.2) 

where a is the inertial moment of the ski with respect to its baricentre, m s its mass and 
m the mass of the skier, so that M = m + m s . In Bressan (1991) it was assumed that: 

uc < 1 -uG[r_,r_|_] (5-3) 

dl 
S 9n(-Q-) = -sgn(c) (5.4) 

where, by definition, sgn(0) = 0. It is easy to verify that (5.3,4) hold under the hypotheses 
(5.2). 

We now consider the case in which the curvature c is constant. Then X s = and (5.1) 
reduces to: 

( It \ i i i 2 1,2 2 (2n+l)(n+l) , (n+l)\ _1 

s = p {ac''(s) + m s + m + mc' : 'b z u z '± ^ '- - m c b u - '- ) /g ^ 

p = —M g (1 — c u) y'(s) 

Now we want to verify that, under suitable conditions on the curvature, the control system 
(5.5) satisfies the hypothesis (PI) of the third section. If c = then the control disappear 
and we obtain a dynamical system. If p = or y'(s) = the set of velocities is a segment 
and we are done. Hence from now on we assume that c,p,y'(s) ^ 0. From the second 
equation of (5.5) we can express u as a function of p and substitute its value in the first 

equation obtaining: 

P 

s 



P{p) 
where: 



. 2 l 'mb 2 (2n + l)(n + l)\ . fmb[b(2n+ 1) - 3n] (n+ 1) 



P(P)=F .^or,,!./ +P 



M 2 g 2 [y'(s)] 2 Qn 2 J \ Mgy'(s)3n 
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6 a c 2 n 2 + m b 2 (2n 2 + 3n + 1) — 6 m b n (n + 1) + 6 (m + m s ) n 2 



Qn 2 
The discriminant of the polynomial P is: 

_ mb 2 (n+l)[2ac 2 (2n + l) + m(n + l) +2m s (2n + l)] 
{ > ~ 3M 2 g 2 [y'(s)] 2 n 2 

then P(p) has no zeros. 

The zeros of the second derivative of s with respect to p are: 

. _ gMy'(s)[3n-b(2n + l)] 
P± ~ b(2n + l) 

± ^^^J — V2(ac 2 + m s )(2n 2 + 3n+l)+m(n 2 -l). 
Vm b (2n + 1) (n + 1) 



Then we compute u± G M in such a way that p(u±) = p±: 



n 
u± = — 



3 v/m (n + 1) ± y/2 (a c 2 + m s ) (2 n 2 + 3 n + 1) + m (n 2 - 1) 



bc v ^m(2n + 1) (n+ 1) 
If c < then «_)_ < and the condition -u_ < gives: 



2 m (4 71 + 5) 777-4 



c 2 < ^1A^^_1L _ — . (5. 6 ) 

a(2n + l) a 



If c > then (5.6) gives U- > but if we assume: 

3n 



r + > (5.7) 



6c(2n + l) 
then the condition U- > r+ gives: 

ci2 c 2 + a\ c + ao > (5.8) 

where: 

a 2 = (2n + l)(n + l) [mb 2 r+ (2n + 1) (n + 1) -2 an 2 ] (5.9) 

a\ = —6mbr+ n (2n + 1) (n+ l) 2 ao = 2 n 2 (n + 1) [m (An + 5) — m s (2n + 1)]. 

For ?7 sufficiently large and for standard data, the conditions (5.7,8,9) are similar to the 
condition (5.3) only slightly more restrictive. If we assume that (5.6,7,8,9) are verified 
then s(ft(u)),u G [r+,r_] is convex or concave, depending on the signes of p, c. Moreover 
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we have that s(p) tends to zero as \p\ tends to infinity and it is easy to check that (PI) is 
verified. 

The assumptions (5.6,7,8,9) are natural in fact in the application we have some bounds 
on the curvature. 

We can now compute the associated system with control appearing linearly, as in 
section 3, to obtain information about the time optimal controls. If c = the control does 
not appear and we have a dynamical system as observed above. Hence from now on we 
assume that c^O. We have: 

M 9 ,,w* +I -, , J G = [ M 9 " +I 7 J (5-10) 

— 2 y («) ( 2 - c (r+ + r_)) / \ ^ y (s) c (r+ - r_) 

where X± = X(r±). After straightforward calculations it follows: 

L J 2 v-ps/'OO 



where: 



X + +X_ X+ — X_ L ± H r ,-r_)- J(2-c(r + + r_)) 

2X + X_ 2X + X_ v + y v v + yy 



and: 



Therefore the equation for turnpikes is: 



A B = -^(-2JLp 2 !/ "( S ) +m5 Lc(^)) 2 (r + -r_)). 



2 _ ^ffc(r + -r_) (y'(s)) 2 

p " 2j ^r (5 - n) 

Observe that the first factor of (5.11) is negative in fact sgn(J) = —sgn(c) (from (5.4)). 
For example if we consider the curve: 



/ cas(\c\s) sin(\c\s) \ g g 

V |q |c| / 



we have: 



l mg(r + — r-) cos 2 (\c\ s) 



2 J sirac s) 



We can argue as in section 4 for time optimal control problems. Given a minimum 
time problem, if every time optimal trajectory for the system (2. 4), (5. 10) contains a 
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Z— trajectory then the time optimal control for the ski does not exist. If the opposite 
happens then the time optimal control exists and is bang-bang. The geometric structure 
of the system (2. 4), (5. 10) is very similar to that one of the system (2. 4), (4. 4). The only 
difference is the sign of the function / of (2.14) that can be different. Consider again the 
curve (5.12). Since we are considering a skier, (5.12) makes sense only if we consider the 
restriction s E [0, n/\c\] if c < and s E [ir/\c\, 2n/\c\] if c > 0. 

For the case c > we can repeat the reasoning made for the swing, that is for the system 
(2. 4), (4. 4). Indeed, with the above restriction, the system (2. 4), (5. 10) is similar to the 
system (2. 4), (4. 4) restricted to the set Q of section 4. In particular we have that the time 
optimal control to reach a given position s, with fixed initial data, in minimum time does 
exist and is bang-bang. The natural hypothesis is that the skier has at least one choice 
that permit him to reach s without turning back, this means s(t) > for every t. Hence 
we have that the optimal trajectory lies in the part of the plane {(s, p) : p > 0}. Therefore 
the time optimal control has at most two switchings and can be determined as in section 
4. 

If c < we are in the opposite case because the natural restriction force us to stay in 
the part of the plane in which there are turnpikes. The time optimal control does not 
necessarily exist. Again we can repeat the reasoning of section 4. 

For Mayer problem (4.9) we have obtained the nonexistence assuming some conditions on 
r±, but it can happen that these conditions are not phisically acceptable for the ski model. 
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